Mon. Not. R. Astron. Soc. 000, 1-14 () Printed 21 June 2011 (MN WF^X. style file v2.2) 



Sunyaev Zel'dovich observations of a statistically complete 
sample of galaxy clusters with OCRA-p 

Katy Lancaster 1 , Mark Birkinshaw 1 , Marcin P. Gawronski 2 , 

Richard Battye 3 , Ian Browne 3 , Richard Davis 3 , Paul Giles 1 , Roman Feiler 2 , 

Andrzej Kus 2 , Bartosz Lew 2 , Stuart Lowe 3 , Ben Maughan 1 , Abdulaziz Mohammad 1 , 

Bogna Pazderska 2 , Eugeniusz Pazderski 2 , Mike Peel 3 , Boud Roukema 2 

and Peter Wilkinson 3 

1 H. H. Wills Physics Laboratory, University of Bristol, Tyndall Avenue, Bristol BS8 1 TL 

2 Toruri Centre for Astronomy, Nicolaus Copernicus University, ul. Gagarina 11, 87-100 Torun, Poland 

3 Jodrell Bank Centre for Astrophysics, The University of Manchester, Alan Turing Building, Manchester M13 9PL 



Received **insert**; Accepted **insert** 



ABSTRACT 

We present 30 GHz Sunyaev Zel'dovich observations of a statistically complete sample 
of galaxy clusters with OCRA-p. The clusters are the 18 most X-ray luminous clusters 
at z > 0.2 in the ROSAT Brightest Cluster Sample. We correct for contaminant 
radio sources via supplementary observations with the Green Bank Telescope, also at 
30 GHz, and remove a cluster that is contaminated by an unresolved X-ray source. 
All 17 remaining clusters have central SZ effects with Comptonisation parameter yo 
exceeding 1.9 x 10~ 4 , and 13 are detected at significance ^ 3er. We use our data to 
examine scalings between yo and X-ray temperature, X-ray luminosity, and the X-ray 
mass proxy Yx, and find good agreement with predictions from self-similar models of 
cluster formation, with an intrinsic scatter in yo of about 25%. We also comment on 
the success of the observations in the face of the contaminant source population, and 
the implications for upcoming cm-wave surveys. 

Key words: cosmology: observations - cosmic microwave background - galax- 
ies: clusters: individual (A1835, ZWCL1953, A689, ZWCL3146, RXJ1532. 9+3021, 
A2390, A2219, RXJ2129.6+0005, A2261, A781, A697, A1763, A68, A520, A267, 
RXJ0439.0+0715, ZWCL7160, A773) - methods: observational 



1 INTRODUCTION 

The thermal Sunyaev Zel'dovich (SZ) effect (Sunyaev & 
Zel'dovich 1972) is a spectral distortion of the Cosmic Mi- 
crowave Background (CMB) radiation due to inverse Comp- 
ton scattering by the hot gas in galaxy clusters. It has long 
been exploited in cosmological and cluster studies in order 
to derive, for example, the Hubble constant (e.g. Hughes 
& Birkinshaw 1998, Mason et al. 2001, Reese et al. 2002, 
Saunders et al. 2003, Bonamente et al. 2006) and the gas 
mass fraction (e.g. Grego et al. 2001, Lancaster et al. 2005, 
LaRoque et al. 2006). Thanks to well developed techniques, 
detections are becoming routine although signal-to-noise 
remains quite poor. However, we are entering an era of 
purpose-built instruments so this is set to improve dramat- 
ically, enabling SZ research to reach its evident potential. 



The main focus of the SZ community at present is to 
utilise the redshift-independence of the SZ surface bright- 
ness in order to perform blind surveys for galaxy clus- 
ters. While other techniques suffer from large intrinsic bi- 
ases and complex selection effects, SZ surveys will pro- 
duce almost mass-limited catalogues and thus far superior 
datasets for constraining cosmological models. The dedi- 
cated SZ surveys, for example Planck (Ade et al. 2011), SPT 
(Staniszewski et al. 2009, Vanderlinde et al. 2010), ACT 
(Menanteau et al. 2010, Marriage et al. 2010) and the SZA 
(Muchovej et al. 2011) are now generating results. Many 
more are expected in the near future, e.g. from AMI (Zwart 
et al. 2008). In order to fully exploit the results of these sur- 
veys, it will be necessary to improve understanding of both 
the 'selection effect' due to the presence of unsubtracted ra- 
dio sources, and also the scalings between cluster SZ observ- 
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ables and various physical quantities, especially the cluster 
mass. Various groups have undertaken similar studies (e.g. 
Benson et al. 2004, Morandi et al. 2007, Bonamente et al. 
2008, Huang et al. 2010). To make further advances, reliable 
observations of large, well-selected samples are required. 

The One Centimetre Receiver Array prototype, OCRA- 
p, is a two-element receiver mounted on the 32-m telescope 
at the Toruh Centre for Astronomy of the Nicolaus Coperni- 
cus University, Poland. OCRA-p proved its SZ capabilities 
by detecting four well known clusters at high significance 
(Lancaster et al. 2007). We present the first in a series of 
papers discussing an X-ray selected sample of 33 clusters 
with well understood selection effects. This paper contains 
details of our OCRA-p observations of a statistically com- 
plete subsample of 18 clusters. The structure of the paper 
is as follows. In Section 2 we give a brief description of the 
Toruh 32-m telescope and the OCRA-p receiver. Section 3 
contains details of the cluster sample. The observations and 
data reduction are described in Section 4, and the problem of 
radio source contamination is detailed in Section 5. Section 
6 gives a brief overview of our X-ray analysis, and Section 7 
presents our results. We discuss these results and conclude 
in Section 8. Throughout the paper we adopt the follow- 
ing cosmological parameter values: Ho — 70kms~ 1 Mpc , 
n m0 = 0.3, n A0 = 0.7. 



2 THE TORUN TELESCOPE AND OCRA-P 

The Toruh observatory is located in Piwnice, 15 km outside 
Toruh in northern Poland. The telescope consists of a 32-m 
parabolic dish and 3.2-m sub-reflector, with a fully steerable 
classical Alt-Az mount. It has receivers operating at 1.4-1.7, 
5, 6.8, 12 and 30 GHz. The telescope is used for a variety 
of studies including interstellar molecular spectroscopy (e.g. 
Blaszkiewicz & Kus 2004) and VLBI (e.g. Bartkiewicz et al. 
2005). 

OCRA is a planned 30 GHz 100-element continuum 
receiver (see Browne et al. 2000). A prototype receiver, 
OCRA-p, funded by a grant from the Royal Society Paul 
Instrument fund, has been on the telescope since 2004. The 
instrument is described in detail in Lowe (2005) and Lowe 
et al. (2007); here we present a short summary. The basic 
radiometer design is based on the prototype for the Planck 
Low Frequency Instrument (LFI, Mandolesi et al. 2000), and 
is similar to the WMAP K-band receivers (Jarosik et al. 
2003). OCRA-p provides only two horn-feeds, the beams of 
which are separated by 3'.1 and have FWHM 1'.2. As the 
beam separation is small, it is possible to reduce the effects of 
atmospheric and gain fluctuations by switching between the 
beams and taking the difference. This can be improved upon 
by further levels of switching. The full switching strategy for 
SZ observations is described in Section 4.1. In addition to 
the SZ program, OCRA-p has been used to study planetary 
nebulae (Pazderska et al. 2009) and for various radio source 
studies including flux density measurements to support the 
Very Small Array CMB program (see Gawrohski et al. 2010 
for details of the OCRA measurements, and Genova-Santos 
et al. (in preparation) for the details of the CMB work) . The 
next phase of the project is an 8-element array, OCRA- 
Faraday (Peel 2009), which was mounted on the telescope 
in late 2009. This is still in the commissioning phase but is 



expected to be fully operational in autumn 2011, and can 
be upgraded to 16 elements in the future. 



3 THE CLUSTER SAMPLE 

In order to derive meaningful cosmological results, or indeed 
to comment on universal cluster properties, it is necessary to 
use a 'fair' sample, or at least one where the selection biases 
are clear and well understood. Large samples fulfilling these 
criteria are rare in SZ astronomy, largely due to the various 
practical challenges which observers must overcome. In this 
work, we attempt to select and observe such a sample, and 
turn to the ROSAT All-Sky Survey in the first instance. 

The Brightest Cluster Sample (BCS, Ebeling et al. 
1998) is derived from a careful analysis of the ROSAT All- 
Sky Survey (RASS, Voges et al. 1999) data at 5 > 0°. 
The resulting catalogue is 90% flux complete to a limit 
of fx = 4.4 x 10~ 12 erg cm -2 s" 1 in the 0.1-2.4keV energy 
band. Since the SZ surface brightness of a cluster is essen- 
tially unaffected by redshift, we select on luminosity, impos- 
ing a limit of Lx ?S 13.08 x 10 44 ergs _1 (the luminosity of 
A773 in the BCS catalogue). In addition, we impose the cri- 
terion z > 0.2 due to the 3'.1 beam-throw of OCRA; obser- 
vations at lower redshift are inefficient with this instrument. 
Our basic sample then consists of 18 clusters. We do not re- 
ject any clusters based on their radio-source environments at 
this stage, even though we appreciate the difficulties which 
this policy may introduce. This has the advantage of allow- 
ing us to investigate any correlation between gas properties 
and those of a central radio galaxy. 

We reject one cluster from the basic 18-cluster sam- 
ple, A689, as on further investigation we have found that 
the ROSAT data are contaminated by a bright central 
point source (Giles et al., in preparation). After removing 
this source from the Chandra data, the luminosity is just 
(1.66 ± 0.24) x 10 44 ergs _1 , more than a factor 8 below the 
limiting luminosity of the cluster sample. The sample is sum- 
marised in Table 1. These clusters appear not to be excep- 
tional in their X-ray structures, and, indeed, follow the ex- 
pected X-ray scaling relations (see Section 7.3). 



4 OCRA DATA 

4.1 Observing Strategy 

The clusters were observed in good weather periods between 
September 2005 and January 2009. We aimed to achieve a 
uniform noise level < 0.5 mjy, which required ~ 220 minutes 
of data at each pointing centre. Trail fields, blank patches of 
sky separated from the target fields by 15-20m in RA, were 
observed over the same range of hour angle after each clus- 
ter observation, using the same position switching routine 
(described below) as for the cluster fields. The integration 
times for the target and trail fields are identical. 

For each observation, whether the cluster or its corre- 
sponding trail field, we employ a position switching strat- 
egy to remove atmospheric contamination from our data. 
The feeds are positioned such that one beam, beam B, is 
coincident with the cluster centre and the other, beam A, 
provides a measure of the blank sky signal. (For extended 
sources, beam A may measure a small signal itself. This 
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Table 1. Initial cluster parameters. Coordinates correspond to the OCRA pointing centres which in some cases were shifted slightly to 
avoid bright point sources. ARA and ADec are the offsets from the maxima of the X-ray surface brightness in the Chandra data. Lx 
values refer to the 0.1-2.4 keV ROSAT energy band and are taken from Ebeling et al. (1998) along with the redshifts. A689 was originally 
in the sample, but on closer inspection it was found to have significant X-ray point source contamination. Once this was removed, the 
cluster luminosity fell to the bracketed value, well below our selection limit of 13.02 X 10 44 ergs -1 . 



Cluster 


RA 


Dec 


ARA 


ADec 


z 


L x 




(J2000) 


(J2000) 


(arcmin) 


(arcmin) 




(lO^ergs" 1 ) 


A1835 


14 01 01.99 


02 53 12.8 


0.02 


0.53 


0.25 


38.53 


ZWCL1953 


08 50 03.00 


36 04 16.0 


-0.80 


-0.04 


0.32 


34.12 


A689 


08 37 25.01 


U 59 40.9 






0.28 


30.41 (1.66) 


ZWCL3146 


10 23 35.98 


04 11 56.0 


-0.90 


0.75 


0.29 


26.47 


RXJ1532.9+3021 


15 32 58.99 


30 21 11.2 


1.13 


0.20 


0.35 


24.40 


A2390 


21 53 34.61 


17 40 10.9 


-0.50 


-1.56 


0.23 


21.44 


A2219 


16 40 21.80 


46 42 47.8 


0.24 


0.30 


0.23 


20.40 


RXJ2129.6+0005 


21 29 40.50 


00 05 49.9 


0.11 


0.51 


0.24 


18.59 


A2261 


17 22 27.00 


32 07 04.1 


-0.03 


-0.88 


0.22 


18.18 


A781 


09 20 25.99 


30 29 57.8 


-0.04 


-0.11 


0.30 


17.22 


A697 


08 42 58.01 


36 21 45.0 


0.03 


-0.20 


0.28 


16.30 


A1763 


13 35 23.21 


41 00 04.0 


0.04 


0.08 


0.22 


14.93 


A68 


00 37 05.30 


09 09 10.8 


-0.20 


-0.38 


0.26 


14.89 


A520 


04 54 19.01 


02 56 48.1 


0.00 


0.00 


0.20 


14.44 


A267 


01 52 42.00 


01 00 25.9 


-0.02 


-0.24 


0.23 


13.71 


RXJ0439.0+0715 


04 39 01.01 


07 16 54.8 


0.08 


0.85 


0.23 


13.25 


ZWCL7160 


14 57 19.99 


22 20 35.9 


1.13 


0.03 


0.26 


13.19 


A773 


09 17 51.00 


51 43 19.9 


-0.32 


-0.34 


0.22 


13.08 



must be properly accounted for - see Section 6.2.) We switch 
between the A and B beams at a rate of 277 kHz, record- 
ing the A — B difference every second. We integrate in this 
position for 25 seconds. The telescope then slews to reposi- 
tion the beams such that beam A measures the cluster and 
beam B measures the sky background, and the differencing 
is repeated, this time integrating for 50 seconds. We then 
return to the starting position and integrate for another 25 
seconds. The differenced signals from the first and third po- 
sitions, {A\ — Bi) and (A3 — B3) are summed, and then sub- 
tracted from the differenced signal from the second position 
(A2 — 82). This recovers twice the cluster signal relative to 
the background regions. The symmetric nature of the switch- 
ing pattern ensures proper subtraction of most atmospheric 
emission and other contaminants such as ground pick-up, 
provided that the target region is not rapidly changing in 
azimuth and elevation. 



4.2 Calibration 

The data are calibrated against an internal noise source, 
which is itself calibrated via observations of the well-known 
bright radio source NGC7027 of flux density 5.37 Jy. We ob- 
serve significant changes in the level of the noise source in 
the form of: (a) sudden 'jumps', sometimes due to the tele- 
scope reconfiguration and (b) additional smaller amplitude 
fluctuations which appear random in nature. We note that 
other than the low-amplitude random fluctuations, the volt- 
age level is generally stable on timescales of a few days. By 
performing regular calibration (~ daily) we reduce the ef- 
fects of the larger fluctuations. In addition, we adopt a con- 
servative approach by comparing the total powers for the 
NGC7027 and cluster observations and rejecting any data 
for which they differ significantly, although a study suggests 
the calibration often stays virtually the same as the total 



power drifts. Residual uncertainties in the calibration of the 
system at the level of 5% may be expected. 

The telescope pointing is calibrated via azimuth and 
elevation scans across a bright source close to each clus- 
ter. Where an acceptable pointing calibration does not exist 
within 60 minutes of a cluster observation, the cluster data 
are rejected. 

4.3 Statistical Data Analysis 

After combination of the second-by-second average data into 
double-differenced measurements of the brightness of the 
sky, and calibration, the data are examined for periods of 
increased noise (which might arise from receiver instabilities 
or bad weather conditions) or individual anomalous points. 
The combination of the data into final averages was per- 
formed including statistical tests for outlier data. The frac- 
tions of the data points rejected by 3a or 5a cuts were small 
in all cases, and no cut-dependent changes larger than a 
small fraction of the error on the mean were seen. The dis- 
tributions of data values in the double-difference data are 
close to Gaussian, with a slight tendency to show elevated 
wings in the distributions: the estimates of the error on the 
mean (Table 5) take account of the full distribution, and not 
the distribution after truncation of outliers. 



5 RADIO SOURCE CONTAMINATION 

Radio source contamination remains a significant problem 
for observations of the Sunyaev Zel'dovich effect. If we know 
the positions and 30 GHz flux densities of all sources near 
our observing locations, we can correct the data and estab- 
lish the true magnitude of the SZ effect for each cluster. 
Ideally, we would consult a survey of the radio sky at a 
frequency close to 30 GHz in order to identify contaminant 
radio sources. The only such survey currently available is 
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the WMAP point source catalogue (Wright et al. 2009) , but 
this is complete only to the 2 Jy level and so is three or- 
ders of magnitude too shallow for our purposes. We consult 
the all-sky NVSS (Condon et al. 1998, 1.4 GHz) and GB6 
(Gregory et al. 1996, 4.85 GHz) catalogues and retain all 
sources within 5' of the pointing centre for each cluster. In 
addition, we include any sources found in similar projects 
and reported in the literature. The list of 58 contaminating 
sources is presented in Table 2. 

5.1 GBT observations 

We observed the 58 sources over the course of several 
sessions at the Green Bank Telescope (GBT) in January 
2008 (observations limited due to poor weather conditions), 
then January and May 2009 (re-scheduled time, again com- 
promised in part by non-ideal conditions). We used the 
GBT Ka-band receiver and the Caltech Continuum Back- 
end (CCB). For a detailed description of this configuration 
of equipment and the observing methods used, see Mason 
et al. (2009); we will describe the salient points. 

5.1.1 Calibration 

Every observing session included at least one flux calibration 
observation, using the most accessible of the GBT standard 
calibrators 3C 147, 3C 48 or 3C 286. We observed sources 
cluster by cluster, and thus they tend to lie in small regions. 
The GBT pointing is stable for up to 40 minutes for night 
time observing, but we rarely required re-pointing near a 
field because slews tended to occur after only a few source in- 
tegrations. We chose one pointing calibrator per cluster (and 
hence its associated sources). We took pointing calibrators 
from the GBT database, choosing the brightest available 
within 10 degrees of the cluster centre. We also performed 
at least one sky dip per observing session to establish the 
sky temperature and opacity. 

5.1.2 Differencing 

We employed the standard GBT nodding strategy, which is 
similar to the OCRA differencing scheme, in order to re- 
move atmospheric contamination. A complete nod (ie the 
sequence of beam 1 on source, beam 2 on source, beam 2 
on source, beam 1 on source) takes 1.5 minutes. Each source 
was observed for at least one complete nod. The number of 
nods was governed by telescope scheduling, though we tried 
to include extra nods for sources expected to have lower flux 
densities. 

5.1.3 GBT flux density determinations 

The CCB has a 14 GHz bandwidth, split into four frequency 
sub-bands centred at 27.75, 31.25, 34.75 and 38.25 GHz. We 
can thus measure simultaneously four separate flux densities 
for each source, from which we can derive a 30 GHz value by 
fitting a power law to the available data. We use a Monte- 
Carlo method to estimate the uncertainties on the interpo- 
lated flux densities. For poor-quality data, we were unable 
to derive accurate 30 GHz flux densities from the GBT data 
alone, as discussed in Section 5.2 



5.2 Determining the 30 GHz flux densities 

The quality of our GBT data is rather mixed due to the wide 
range of observing conditions experienced. As most sources 
were only observed once, we are also unable to constrain 
source variability with GBT data alone. In order to supple- 
ment our 30 GHz measurements, we make use of the flux 
densities reported in Coble et al. (2007) who observed all of 
our clusters except Abell 2390 at 28.5 GHz with the OVRO 
and BIMA interferometers. Such interferometric data not 
only measure the SZ effect but also enable identification of 
point sources via the longest baselines, where the SZ signal 
will be negligible. The uncertainties on the Coble et al. flux 
densities are generally a few tenths of a mjy, so given the 
proximity of the two observing frequencies, we are able to 
place tighter constraints on our 30 GHz source flux densities 
by fitting to both datasets, which is particularly useful for 
sources with poor GBT data. In addition, we are able to 
identify variable sources and estimate the effect on our SZ 
measurements. 

Most source flux densities are well constrained by a 
combination of GBT and OVRO/BIMA data. We fit a power 
law to the measurements in the four GBT frequency sub- 
bands, plus the interferometer data at 28.5 GHz. For sources 
with good GBT data, the Coble et al. points serve as a valu- 
able check for variability, but otherwise have minimal effect. 
In the three cases where significant variability is found, we 
take the two differing flux density measurements, and use 
half the difference as an estimate of the uncertainty in our 
measurement due to variability. We then add this in quadra- 
ture to the measurement error. The source exhibiting the 
greatest variability, source 1 in Abell 2261, has little effect 
on our results since it lies well away from the cluster centre 
and data where it lies close to a reference arc are flagged out 
(Section 5.3). Where no 28.5 GHz flux density is available, 
we proceed with the GBT data alone. Such cases are un- 
likely to be crucial, as the non detection by OVRO/BIMA 
suggests the source is either rather faint, or well away from 
the central regions of the cluster. Many sources have GBT 
data which are usable but rather noisy. In such cases, the 
28.5 GHz data are invaluable in tying down the 30 GHz flux 
density due to their comparatively small uncertainties. 

Seven sources have no usable GBT data: two in A1835, 
one in A2219, one in A1763, one in A68 and two in A773. 
We now discuss our method for each source in turn. Source 
1 in A1835 was measured twice by OVRO/BIMA, and also 
appears in NVSS. We fit a power law to the three mea- 
surements. As our observing frequency is close to that of 
OVRO/BIMA, assuming a constant spectral index is un- 
likely to have a large impact. We rule out significant variabil- 
ity. Source 2 in A1835 was also measured twice at 28.5 GHz, 
but does not appear in NVSS. We consult FIRST for a 
lower frequency measurement, and proceed as for source 
1. For source 2 in A2219, we again fit a power law to the 
OVRO/BIMA and NVSS datapoints. Source 6 in A68 does 
not appear in any of the lower frequency catalogues but is 
well measured by OVRO/BIMA. It is outside the OCRA 
beams so its flux density has only a small effect on the SZ 
measurement. Both sources in A773 are faint in NVSS and 
are not detected by OVRO/BIMA. The noise on the OVRO 
map is 0.078 mjy, giving a 5a upper limit of 0.39 mjy. As 
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Table 2. Radio sources found within 5 arcmin of the cluster pointing centre in the NVSS (Condon et al. 1998) and GB6 (Gregory et al. 
1996) surveys, and also in the observations of Coble et al. (2007). Notes: All sources were observed with the GBT unless marked (*). 
30 GHz flux densities were derived from our GBT measurements plus supplementary information at other frequencies as indicated in 
the final column (OB for OVRO/BIMA, N for NVSS, G for GB6 and F for FIRST). Quoted errors arc ±la (68% confidence interval). 
Some sources were observed more than once with OVRO/BIMA. Sources showing evidence for variability are marked (V). Sources with 
apparently negative 30 GHz flux densities are shown in italics: these data were not used to correct the SZ measurements. 



Cluster 




Src 


RA 


Dec 


Si. 4 


S4.85 


S28.5 


S30 


Notes 






# 


(J2000) 


(J2000) 


(mJy) 


(mJy) 




(mJy) 




A1835 




1* 


14 01 02.1 


02 52 41.0 


39.3±1.6 


- 


3.31±0.14 
2.88±0.07 


2.93±0.18 


OB N 






2* 


14 01 00.5 


02 51 53.0 


1.6±0.1 


- 


1.26±0.14 
1.36±0.08 


1.33±0.07 


OB F 


ZWCL1953 




1 


08 50 07.8 


36 04 21.7 


19.4±1.2 


- 


- 


0.79±0.70 


GBT 






2 


08 50 13.4 


36 04 22.4 


20.8±1.0 


19.0±4.0 


1.19±0.17 


1.28±0.16 


GBT OB 


ZW3146 




1 


10 23 37.2 


04 09 06.4 


31.5±1.0 


- 


2.03±0.22 
2.12±0.15 


2.35±0.81 


GBT OB 






2 


10 23 39.6 


04 11 15.4 


7.1±0.5 


- 


0.41±0.07 


0.33±0.07 


GBT OB 






3 


10 23 45.1 


04 10 40.7 


95.8±3.4 


42.0±7.0 


5.35±0.18 
5.70±0.10 


5.27±0.11 


GBT OB 






4 


10 23 45.2 


04 11 39.7 


3.6±0.4 


- 


0.85±0.10 


0.80±0.10 


GBT OB 


RXJ1532^ 


-0005 


1 


15 32 47.4 


30 18 46.0 


18.0±1.0 


- 


- 


1.18±0.21 


GBT 






2 


15 32 50.7 


30 19 47.6 


7.9±0.5 


- 


6.58±0.20 


6.56±0.21 


OB N 






3 


15 32 53.8 


30 20 59.8 


22.8±0.8 


20.0±4.0 


3.25±0.18 


3.27±0.12 


GBT OB 






4 


15 32 54.31 


30 23 01.5 


4.4±0.4 






0.92±0.21 


GBT 






5 


15 33 03.42 


30 23 47.2 


2.6±0.4 


- 


- 


1.12±0.50 


GBT 


A2390 




1 


21 53 32.4 


17 42 19.8 


20.9±1.6 


- 


Not observed 


0.23±0.12 


GBT 






2 


21 53 36.8 


17 41 44.8 


235.3±8.3 


220±20 


Not observed 


45.77±0.10 


GBT 






3 


21 53 40.3 


17 42 56.7 


12.2 ±1.0 


- 


Not observed 


-0.14±0.09 


GBT 


A2219 




1 


16 40 21.83 


46 42 47.8 


239.1±8.3 


84.8±8.0 


14.87±0.17 


13.74±1.00 


GBT OB 






2* 


16 40 23.83 


46 41 47.3 


7.9±1.0 


- 


0.97±0.17 


0.94±0.17 


OB N 






3 


16 40 14.98 


46 42 28.7 


6.1±0.5 


- 


- 


0.42±0.09 


GBT 






4 


16 39 58.07 


46 40 37.2 


14.1±0.5 


- 


- 


0.96±0.09 


GBT 


RXJ2129.6-I 


-0005 


1 


21 29 36.61 


00 02 35.4 


2.2±0.4 


- 


- 


-0.13±0.21 


GBT 






2 


21 29 40.00 


00 05 22.9 


25.4±1.2 


- 


2.33±0.10 


2.04±0.06 


GBT OB 






3 


21 29 40.14 


00 01 44.4 


4.5±0.5 


- 


- 


0.11±0.13 


GBT 






4 


21 29 55.24 


00 07 56.9 


34.3±1.8 


- 


- 


2.82±0.11 


GBT 


A2261 




1 


17 22 16.9 


32 09 10.5 


23.0±1.5 




9.32±0.22 
10.48±0.16 


15.58±3.89 


GBT OB (V) 






2 


17 22 27.7 


32 07 57.8 


5.3±0.5 






0.13±0.08 


GBT 






3 


17 22 45.3 


32 09 27.0 


4.85±0.5 






0.14±0.10 


GBT 






4 


17 22 06.55 


32 07 01.7 


5.8±0.4 






0.14±0.08 


GBT 


A781 




1 


09 20 08.51 


30 32 14.3 


19.9±0.9 






1.42±0.60 


GBT 






2 


09 20 14.12 


30 29 02.8 


17.9±0.7 






1.22±0.53 


GBT 






3 


09 20 21.75 


30 32 27.0 


2.7±0.5 






0.11±0.18 


GBT 






4 


09 20 22.90 


30 29 45.6 


73.1±2.6 


32.0±5.0 


5.33±0.18 


5.35±0.30 


GBT OB 






5 


09 20 30.83 


30 28 02.1 


15.8±1.6 






-0.07±0.16 


GBT 






6 


09 20 47.35 


30 28 21.3 


5.4±0.6 






0.08±0.17 


GBT 


A697 




1 


08 42 40.22 


36 19 16.4 


5.4±0.5 






1.64±0.54 


GBT 






2 


08 42 59.67 


36 17 43.7 


32.5±1.4 






1.71±0.65 


GBT 
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Table 2 - continued 



Cluster 


Src 


RA 


DEC 


Si. 4 


<S4.85 


S28.5 


S30 


Notes 




# 


(J2000) 


(J2000) 


(mJy) 
(NVSS) 


(mJy) 
(GB6) 


(mJy) 
(O/B) 


(mJy) 
(GBT+) 




A1763 


1 


iS 35 15.60 


41 00 25.8 


7.9±0.5 






-0.10±0.14 


GBT 




2* 


13 35 19.86 


41 00 04.0 


859.2±29.5 


225±20 


31.30±0.41 


29.49±0.38 


OB G N 


A68 


1 


00 36 52.94 


09 05 21.7 


22.8±0.8 


- 


- 


1.57±0.05 


GBT 




2 


00 37 05.21 


09 13 33.6 


4.4±0.4 






0.60±0.05 


GBT 




3 


00 37 06.35 


09 07 30.4 


40.2±1.3 




1.20±0.12 


1.21±0.06 


GBT OB 




4 


00 37 07.71 


09 08 24.0 


59.1±2.2 


- 


1.60±0.10 


1.65±0.05 


GBT OB 




5 


00 37 17.93 


09 06 30.5 


3.2±0.6 






-0.05±0.04 


GBT 




6* 


00 37 07.00 


09 07 58.7 


- 


- 


1.38±0.11 


- 


OB 


A520 


1 


04 54 01.11 


02 57 45.6 


6.3±0.5 


- 


7.83±0.25 
4.42±0.88 


5.75±1.74 


GBT OB (V) 




2 


04 54 03.72 


02 56 01.9 


7.7±1.7 


- 


- 


0.49±0.15 


GBT 




3 


04 54 12.41 


02 57 52.3 


7.1±0.5 


- 


- 


-0.15±0.16 


GBT 




4 


04 54 17.14 


02 55 33.5 


14.1±1.0 


- 


0.84±0.11 
1.09±0.09 


0.90±0.08 


GBT OB 




5 


04 54 21.16 


02 55 01.4 


26.5±1.6 


- 


1.00±0.13 
0.74±0.14 


0.63±0.10 


GBT OB 


A267 


1 


01 52 29.46 


00 59 31.8 


30.0±1.0 


- 


2.75±0.20 


3.12±0.23 


GBT OB 




2 


01 52 54.58 


01 02 08.2 


4.2±0.5 




7.55±0.24 
5.53±0.73 


7.05±1.12 


GBT OB(V) 




3 


01 52 34.35 


01 01 20.5 








0.59±0.32 


GBT 


RXJ0439.0+0715 


1 


04 39 01.26 


07 15 42.6 


30.6±1.4 




1.18±0.16 


0.89±0.22 


GBT OB 


ZWCL7160 


1 


14 57 08.06 


22 20 11.2 


13.2±0.6 




0.95±0.05 


0.90±0.06 


GBT OB 




2 


14 57 08.06 


22 20 11.2 


3.9±0.4 




0.99±0.07 


0.91±0.08 


GBT OB 




3 


14 57 15.18 


22 20 36.0 


16.5±1.3 




0.96±0.04 


0.90±0.06 


GBT OB 


A773 


1* 


09 17 45.39 


51 43 11.2 


2.7±0.5 




< 0.39 




Negligible 




2* 


09 18 01.81 


51 44 11.4 


3.1±0.5 




< 0.39 




Negligible 



the sources in A773 lie well away from the OCRA beams, 
we are confident that their effects on our data are negligible. 



5.3 Source corrections 

Due to OCRA-p's position switching strategy, contaminant 
radio sources can affect the data by producing a positive 
signal when they lie close to the cluster centre, or a rogue 
negative signal when they lie in the reference beam of the 
telescope. We simulate the effects of the sources based on 
their measured positions and flux densities, and the OCRA 
beam response. We then use the resulting effective flux den- 
sity to correct the SZ data. Contamination through sources 
at 30 GHz is generally low for bright clusters, producing only 
small corrections to the SZ data, and thus introduces little 
bias as a result of the limitations of the method employed. 
Although we are confident that residual biases are below the 
level of our errors, we recognise that the absence of 30 GHz 
data is a limitation in some cases, and there may be addi- 
tional issues to do with source variability. For clusters where 



radio sources fall in the OCRA reference beams, we are able 
to perform an additional check by flagging out contaminated 
data based on the range of parallactic angles affected by the 
source in question. We compare the flagged data with the 
corrected data, and adjust our final measurement error ac- 
cordingly. 

We recognise that our strategy may miss sources which 
are important at 30 GHz. The flux limit of NVSS is ~ 
2.5 mJy. If we consider a source of flux density ~ 2.4 mJy 
at 1.4 GHz, and assume a typical spectral index of —0.7, ex- 
trapolating to 30 GHz gives a flux density of 0.34 mJy. Even 
if such a source were to be located in the central regions 
of a cluster, the effect on our measurement would still be 
small relative to random errors. Not all sources have typ- 
ical spectral indices, but we note that in the 16 clusters 
additionally observed by Coble et al. (2007), only 1 'new' 
source, i.e. a source not previously identified in NVSS, was 
found. However, the quality of our knowledge of the popula- 
tion of sources with mJy flux densities at 30 GHz is poor. We 
were able to run an additional check for 'new' sources lying 



7 



in OCRA-p's reference arcs. We binned the data by par- 
allactic angle to look for contaminating sources that might 
produce false SZ effects. Although this subdivides the data, 
and so leads to a noise level that is higher by a factor of 
2 — 3 than the central measured flux density averaged over 
all parallactic angles, no previously uncatalogued sources 
were detected. 



6 X-RAY DATA 

6.1 Deriving cluster parameters 

Our sample has 15 clusters in common with an upcom- 
ing X-ray paper, Giles et al. (in prep), some results from 
which were made available in advance for this work. The 
remaining two clusters (ZWCL1953 and RXJ1532. 9+3021, 
which lie outside Giles et al.'s specified redshift range of 
0.15 < z < 0.30) were analysed using an identical method. 
The authors use Chandra ACIS-I and ACIS-S imaging and 
spectroscopy for each cluster as available via the archive. For 
the data preparation, they follow closely the method out- 
lined in Maughan et al. (2008). The Chandra analysis pack- 
age used was CIAO 1 version 4.2 and calibration database 
CALDB 4.1.4. As in Maughan et al. (2008), blank sky, rather 
than local, backgrounds are used because the clusters fill a 
large fraction of the field of view. 

We fit a gas density model by converting the observed 
surface brightness profile into a projected emissivity pro- 
file, which was then modelled by projecting the modifed /3- 
model of Vikhlinin et al. 2006 along the line of sight (see 
Maughan et al. 2008 for details). Gas parameters and er- 
rors were determined from Monte Carlo realisations of the 
projected emissivity profile. At each data point a new ran- 
domised point was drawn from a Gaussian distribution cen- 
tered on the model value at that point, with a standard 
deviation equal to the fractional measurement error on the 
original data point, multiplied by the model value. The er- 
ror bar on this new randomised point was set to the same 
fractional size as the error bar on the original observed point. 

The cluster temperature, gas mass and Rsoo (the radius 
enclosing a mean density of 500 times the critical density at 
the cluster's redshift) were then determined iteratively. The 
procedure followed was to extract a spectrum from within 
an estimated -R500 (with the central 15% of that radius ex- 
cluded), integrate the gas density profile to determine the 
gas mass within the estimated -R500, and thus calculate Yx 
which is the product of kT and the gas mass, and a low scat- 
ter proxy for total mass (Kravtsov et al. 2006). A new value 
of -R500 was then estimated from the Yx-M scaling relation 
of Vikhlinin et al. (2009), 



M 5 



E(z) 



- 2 / 5 A 

-Aym 



3 x 10 14 A/ Q keV 



(1) 



with Aym = 5.77 x 10 14 /i 1/2 Af© and Bym = 0.57. Here, 
M500 is the mass within -R500 (allowing -R500 to be trivially 
computed), and 



E(z) = ^/Q M (1 + z) 3 + (l - fi M - 00(1 + z) 2 + Oa (2) 

describes the redshift evolution of the Hubble parame- 
ter. This Yx~M relation assumes self-similar evolution (as 
E~ 2 ^ 5 ), which has been shown to be a good description of 
observed clusters to z « 0.6 (Maughan 2007). The process 
was repeated until -R500 converged. In these spectral fits, 
the source emission was modelled with an absorbed thermal 
plasma APEC model in the 0.6-9 keV band, with the ab- 
sorbing column fixed at the Galactic value. Nu values were 
taken from the HI map by Kalberla et al. (2005). 

Our final cluster parameters, with Tx, Lx and Yx de- 
termined within a radius of -R500 are given in Table 3. 



6.2 /3-model fitting 

As mentioned in Section 4.1, the small separation of the 
OCRA-p beams relative to the extent of a cluster's SZ decre- 
ment implies that we subtract out a fraction of the SZ signal 
due to our differencing method. This significantly reduces 
the measured signal compared to the intrinsic central sig- 
nal. This is further exaggerated in cases where we have ad- 
justed the pointing centre to avoid bright radio sources. In 
order to remove this effect and recover the 'true' SZ signal, 
we model the cluster atmosphere and estimate the fraction 
of the signal 'missed'. (Note: this analysis was conducted 
purely to correct the SZ measurements and is separate from 
the derivation of X-ray cluster parameters discussed in Sec- 
tion 6.1). In our previous work (Lancaster et al. 2007), we as- 
sumed that the gas density distribution, n e (r) was well fit by 
a standard spherical /3-model (Cavaliere & Fusco-Femiano 
1976): 



,{r) — n c0 I 1 + 



-3/3/2 



(3) 



Here, we improve on this approach by also considering the 
double /3-model, which may provide a better fit to clusters 
which exhibit peaked central emission due to their cool cores. 
We follow Bonamente et al. (2006) in adopting the form 



n e (r) — n a o 



V 1 + 



-301/2 



+ (1-P) 1 + 



-30 2 /2' 



(4) 

although in our analysis, we allow different /3 for the two 
components. Using the SHERPA 2 modelling and fitting ap- 
plication for each set of cluster data, we fit the 1-D radial 
profiles of the X-ray surface brightness by three models: 

(i) The standard single /3-model, (3) 

(ii) A double /3-model, (4), with j3i = /3 2 

(iii) A double /3-model, (4) , allowing /3i and P2 to assume 
different values. 

Chi squared (y 3 ) statistics are used to compare the fits, since 
each radial bin is based on a large number of X-ray counts. 
For each cluster, we accept the simplest model for which 
the fit statistic implies a probability of getting a greater 



1 http://cxc.harvard.edu/ciao/ 



2 http://cxc.harvard.edu/shcrpa/ 
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Table 3. Chandra X-ray parameters for the sample, derived inside -R500 with the central 15% excluded. Lx is the bolometric luminosity. 
Quoted errors are ±1<t (68% confidence interval). 



Cluster 


T x 


L X 


Y x 


^500 




(keV) 


(lO^crgs" 1 ) 


(lO 15 M keV) 


Mpc 


A1835 


4 i+0.32 


2.14±0.01 


1 20+ ' 05 
1 - zu -0.04 


1.40 


ZW1953 


fi nK+0-51 
°- UO -0.50 


1.45±0.05 


U - 58 -0.05 


1.19 


ZW3146 


' - ua -0.37 


1.62±0.02 


77 +0.04 
u -' '-0.04 


1.27 


RXJ1532. 9+3021 


5 62+ - 72 
o,DZ -0.58 


1.23±0.05 


o 45+ 006 

u - 4o -0.06 


1.11 


A2390 


W.20t° f 5 


2.87±0.05 


1 7Q+OO6 
± -' 3 -0.07 


1.53 


A2219 


10.33l°« 


3.32±0.04 


2 m+ 09 

z - ±u -0.09 


1.58 


RXJ2129.6+0005 


fi n 4 +0.53 


1.03±0.04 


n ko+0.05 
u ' oz -0.06 


1.20 


A2261 


6 56+ - 28 
D.OD_q 2g 


1.46±0.02 


n --+0.03 
u -' '-0.03 


1.31 


A781 


t- 09+O.66 
-0.43 


1.02±0.05 


n rcq + 0.06 
u - o,3 -0.04 


1.18 


A697 


Q 99+O.66 
M - z -0.65 


2.70±0.05 


1 ==+0.11 

• L * oo -0.11 


1.45 


A1763 


7 O-1+0.47 
' ■ ol -0A7 


1.48±0.03 




1.34 


A68 


7 14+O.88 
1 ilyi -0.78 


1.09±0.04 


60+ 07 

u ' DU -0.07 


1.23 


A520 


f, 7 o+0.29 


1.49±0.02 


n 74+0 03 

u -' 4 -0.02 


1.31 


A267 


4 53+ - 49 
^■ JO -0.38 


0.75±0.05 


29+ 03 
u - z3 -0.03 


1.08 


RXJ0439.0+0715 


t- 09+O.38 


1.58±0.04 


0.57iH 4 4 


1.23 


ZW7160 


4 64+°- 18 


0.69±0.01 


0.28+Toi 


1.06 


A773 


fi 7^+ - 40 
,o -0.26 


1.13±0.02 


o.65tg;° I 


1.27 



X 2 by chance, Q ^ 6 per cent. For some clusters, the gas 
distribution is sufficiently disturbed that Q < 6 per cent for 
all models fitted: in such cases we adopt the best fit. This is 
not ideal, but we expect the additional error to be at a low 
level in comparison with the random errors on the SZ data. 
We do not observe any systematic effect on our results due 
to the type of /3-model used to correct the SZ data. 

/3-model parameters are presented in Table 4. We re- 
turn briefly to the analysis of Bonamente et al. (2006), 
with whom we have eight clusters in common, for com- 
parison of results. For A697, A68, A267 and A773, both 
groups find that a single /3-model provides an acceptable 
fit, and the derived parameters are in broad agreement. For 
their double /3-model fits, Bonamente et al. fix the two val- 
ues to /3 to be equal, whereas ours were allowed to vary. 
For RXJ2129.6+0005 and A2261, we find that a model 
with /3i = P2 gives the best fit, and the model parame- 
ters are again in broad agreement with Bonamente et al. 
For ZWCL3146, our two /3-values differ but are very simi- 
lar, meaning the overall model fit is in broad agreement with 
Bonamente et al., and indeed our /3i = /3 2 fit is also consis- 
tent. The exception is A1835, where we find the best fitting 
model to be a double /3-model with /3i = 0.675, P2 = 1.059, 

81 = 35.6 arcsec and 62 = 13.2 arcsec, whereas Bonamente 
et al. find /3 = 0.797, 81 = 63.6 arcsec and 82 = 13.2 arcsec. 
Our model with /3i = 02 is also slightly discrepant, as our 
best fit parameters are /3 = 0.695, 9i — 92.1 arcsec and 

82 = 19.1 arcsec. The source of discrepancy may well be in 
the use of differing data sets. We use ACIS-I observations 
totalling 200ks, whereas Bonamente et al. only had access to 
30ks of ACIS-S data which has a significantly smaller field 
of view. 

We use the model gas distribution to simulate an 
OCRA observation. Many clusters were observed at posi- 



tions slightly shifted from their X-ray maxima in order to 
avoid bright radio sources as far as possible. We consult 
our OCRA data for the range of parallactic angles sampled, 
and include the convolution with the OCRA beam response. 
We can thus estimate the fraction F of the 'true' SZ signal 
measured in the real observations, and 'correct' our measure- 
ment (multiply by 1/F) to recover the SZ surface brightness 
that would be measured by an ideal instrument with nar- 
row beams at infinite separation. We note for clusters where 
the double /3-model provides an improved fit, the value of 
F varies from the simpler single /3-model by ~ 10 per cent. 
This change is thus significant for the scientific interpreta- 
tion of our results. Correction factors are presented in Table 
5. Asymmetries in the SZ signal are assumed small in our 
analysis. Checks for such asymmetries by binning the SZ 
data by parallactic angle provide only a low signal-to-noise 
control against such effects which would cause our estimates 
of F to be incorrect. The X-ray isophotes of these clusters 
are not always circular, but are most sensitive to asymmetry 
towards the cluster centre, where a relatively small fraction 
of the SZ effect originates. We estimate that the value of F 
could be systematically uncertain by 20% from this effect. 
This uncertainty dominates over that from varying the den- 
sity model from /3 or double-/3 form to the Vikhlinin et al. 
(2006) form. 

7 RESULTS AND SCALING RELATIONS 
7.1 Sunyaev Zel'dovich Results 

We present our measurements for the Sunyaev Zel'dovich ef- 
fects in our 17-cluster sample in Table 5. For each cluster, we 
first give the raw OCRA flux density as measured directly 
by the telescope, followed by the revised figure once all ra- 
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Table 4. /3-model parameters for the cluster sample. The best fitting model is shown for each cluster, where the parameter p indicates 
the fraction of each component as described by Equation 4. Quoted errors are ±lcr (68% confidence interval). In cases where a double-/3 
model provided the best fit, we choose the larger of the core radii to be r c i The resulting correction factors applied to the SZ data are 
given in Table 5. 



Cluster fix r cl /3 2 r c2 p 

(arcsec) (arcsec) 



A1835 





.675+°' 

— U 


001 

001 


35. 


g+0.6 
—0.6 


1 


0591° 


010 
010 


13 


9+0.1 
— 0.1 





10 


ZWCL1953 





,972^ 


147 
122 


171. 


fi +22.3 
— 21 9 





,9721° 


147 
122 


13 


q+5.7 
— 5 1 


0. 


04 


ZWCL3146 





.7001°, 


008 
007 


23 


g+0.9 
—0.8 





6521° 


032 
(123 


4 


o+0.3 
O no 
— (J.Z 





14 


RXJ1532.9+3021 





.650+°, 


011 
010 


12 







6501°, 


011 
010 


6 


3 +1.6 
— 1.8 





40 


A2390 





■8311° 


010 
009 


92 


Q+1.6 

a -1.6 





8311° 


010 
000 


18. 


c+0.3 
°-0.3 


0. 


10 


A2219 





7251° 


015 
013 


85 


g+6.2 
D -5.0 





7251° 


015 
013 


12 


Q +5.4 
U -5.6 


0. 


58 


RXJ2129.6+0005 





609lo 


011 
011 


21 


4+2-9 
4 -2.5 





609l°j 


011 
011 


5 


n +o.8 


0. 


12 


A2261 





5981° 


003 
003 


30 


7+1.4 
'-1.4 





5981° 


003 
003 


11 


S -0.8 





33 


A781 





.7061°, 


041 
035 


79 


t;+7.5 
°-6.6 
















A697 





6481° 


008 
008 


48 


7+1.5 
' -1.5 
















A1763 





.5861° 


007 
(100 


48 


n +1.6 
u -1.6 
















A68 





7421°' 


003 
002 


51 


c;+3.5 
°-3.3 
















A520 





.8311° 


003 
003 


117. 


rc + 0.6 
°-0.6 
















A267 





6331°, 


(111 
(111 


31 


4+1-5 
1.4 
















RXJ0439.0+0715 





.7351° 


022 
020 


52 


g+4.3 
D -3.9 





7351° 


022 
020 


13 


+1.6 
■ u —1.5 


0. 


23 


ZWCL7160 





.8501°' 


112 
077 


22 


n+4.3 
z -3.4 





.5801° 


005 
006 


7 


9+0.3 
z -0.4 


0. 


09 


A773 





7691° 


041 
035 


129. 


g+14.3 
D -14.4 





7691° 


041 
035 


48 


9+2.7 
-2.6 
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dio source contamination has been corrected for. We then 
give our estimate of the fraction of total signal measured 
due to the instrumental response. Finally, we give the fully- 
corrected temperature decrement expressed in brightness 
temperature units, and the corresponding central Compton- 
isation parameter, yo = —0.19 x (AT/K). We detect 13 
clusters at significance ^ 3<r. 

Our trail fields serve as a valuable test of our observing 
strategy. We present differential and cumulative distribu- 
tions for both the source corrected cluster target and trail 
fields in Figure 1. The flux densities of the trail fields are 
distributed around zero with small scatter, and their distri- 
bution is clearly distinct from that of the SZ data, confirming 
that our observing strategy is effective. 

7.2 SZ— X-ray scaling relations 

According to the self-similar model of galaxy cluster for- 
mation, which assumes evolution solely under gravity, we 
expect simple scalings between various cluster parameters. 
Following Morandi et al. (2007) we derive the following re- 
lations, each of which may be tested using OCRA data. We 
note that while it is advantageous to use the integrated 
Compton parameter for such studies (e.g. Benson et al. 
2004, Bonamente et al. 2008), deriving this from OCRA data 
would introduce further uncertainty since we have measured 
the cluster structures from the X-ray data alone. Thus, we 
restrict ourselves to discussing the central Comptonisation, 
yo- The scalings of interest are first, the relation between 
central Comptonisation and X-ray temperature 

yo (x T* /2 E(z), (5) 



second, that between the Comptonisation and X-ray lumi- 
nosity 

yo cx L^ 4 E(z) 1/4 , (6) 

and finally, that between Comptonisation and Y%, the prod- 
uct of gas mass and X-ray temperature 

yo (X Y^E{zf'\ (7) 

where E(z) is as given in Equation 2. Deviation from the 
expected scalings would suggest contributions from non- 
gravitational processes, and redshift dependence in the slope 
or amplitude of the scalings should provide evidence on the 
nature of these processes. 

We present our scaling relations in figures 2(a) to 2(c). 
For each plot, we colour-code the datapoints based on the 
reliability of the OCRA data. The most reliable points, i.e. 
those with no obvious issues, are shown as red circles. Points 
which have a signal to noise ratio of less than 3 are black 
squares. This is an arbitrary cut-off, and we see no reason 
to exclude these points from our analysis; however as an aid 
to the reader we choose to mark them clearly. Blue crosses 
represent clusters which may suffer from residual source 
contamination, i.e. those which contain a bright (> 3mJy) 
source, for which we have corrected, in their central regions. 
Due to the unavoidable 10 per cent uncertainty in the cal- 
ibration of the GBT data, we cannot be confident that the 
effects of these sources have been accurately removed, and 
thus choose to ignore these three clusters in our subsequent 
analyses. 

The fitting was performed in the linear space of observ- 
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Flux density/ mJy 

Figure 1. Lower panel: histogram of the source-corrected flux densities for the target cluster fields (black solid lines) and the flux 
densities for the reference trail fields (red dashed lines with cross-hatching). Upper panel: the corresponding cumulative distributions for 
the clusters (black diamonds) and the trail fields (red circles). Note that the trail field flux densities arc consistent with zero, suggesting 
that our switching scheme removes atmospheric and ground-based contamination effectively. The cluster and trail field flux density 
distributions clearly differ. 



Table 5. OCRA SZ measurements. 5R aw gives the raw OCRA measurement, whereas Ssc is the SZ flux after correction for contaminant 
point sources. The fourth column lists the measured fractions of the expected SZ signal due to the OCRA switching pattern. AT is 
the brightness temperature decrement, corrected for point sources and the instrumental response. All errors are ztlcr (68% confidence 
interval) . 



Cluster 






•Ssc 


Fraction 


AT 


yo 




(mjy) 




(mJy) 


Measured 


(UK) 


xlO 4 


Af835 


-2.05ztO.46 


-3 


82zfc0.47 


0.363 


-4050zt500 


7.70zfc0.95 


ZW1953 


-1.65±0.37 


-1 


65zfc0.37 


0.404 


-1570zfc350 


2.98zfc0.67 


ZW3146 


-1.06±0.39 


-0 


78zfc0.41 


0.181 


-1660zfc870 


3.15zfcl.65 


RX1532.9+3021 


-1.41zfc0.33 


-1 


llzfc0.34 


0.166 


-2540zfc790 


4.83zfcl.50 


A2390 


-0.78±0.20 


-1 


02zfc0.22 


0.176 


-2230zt480 


4.24zfc0.91 


A2219 


11.10±0.38 


-2 


76zfc0.40 


0.379 


-2810zfc400 


5.34zfc0.76 


RXJ2129.6+0005 


-0.13zfc0.37 


-1 


46zfc0.37 


0.329 


-1710zfc430 


3.25zfc0.82 


A2261 


-3.19ztO.30 


-1 


35zfc0.32 


0.321 


-1610zfc390 


3.06zfc0.74 


A781 


0.68zfc0.42 


-1 


19zfc0.53 


0.330 


-1390zfc620 


2.64zfcl.l8 


A697 


-3.01ztO.46 


-3 


01zfc0.46 


0.422 


-2740zfc420 


5.21zfc0.80 


A1763 


10.58zfc0.80 


-3 


llzfc0.82 


0.296 


-4040zfcl070 


7.68zfc2.03 


A68 


-0.89zfc0.48 


-1 


22zfc0.48 


0.451 


-1040zfc410 


1.98zfc0.78 


A520 


-1.97zfc0.47 


-1 


80zfc0.48 


0.268 


-2580zt680 


4.90zfcl.29 


A267 


-2.73zfc0.43 


-1 


99zfc0.57 


0.471 


-1620zfc470 


3.08zfc0.89 


RX0439.0+0715 


-1.22zfc0.43 


-1 


27zfc0.43 


0.351 


-1390zfc470 


2.64zfc0.89 


ZWCL7160 


-I.10zfc0.34 


-1 


14zfc0.34 


0.188 


-2330zfc700 


4.43zfcl.33 


A773 


-1.20zfc0.32 


-1 


21zfc0.32 


0.382 


-1210zt330 


2.30zfc0.63 
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Figure 2. Scaling relations between the central Comptonisation as derived from OCRA data and various X-ray parameters. We fit 
to the 'robust' and 'low significance' data (red circles and black squares respectively), but omit the clusters with potential residual 
source contamination (blue crosses); the best-fit result is shown as the solid red line. The dotted black line depicts the expectation from 
self-similar evolution. 
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Table 6. Best-fit slopes for the correlations examined in Sections 7.2 and 7.3, along with the expectations from self-similar models for 
comparison. For the SZ/X-ray scalings, we give results for the analysis both with and without corrections for contaminant radio sources 
in the OCRA data. We note that the fits are significantly poorer where the source corrections are omitted. All fits are consistent with 
predictions from self-similar models of cluster formation. 



Correlation 


Expected slope 


Data used 




Best-fit slope 


Fit statistic 


2/0 


1.5 


SZ, X-ray 




1.82 ±0.72 


2.39 






SZ (no source 


correction), X-ray 


1.94 ± 1.10 


6.64 




0.75 


SZ, X-ray 




0.77 ±0.21 


2.46 






SZ (no source 


correction), X-ray 


0.71 ±0.19 


6.86 


yo,Y x 


0.6 


SZ, X-ray 




0.65 ±0.21 


2.50 






SZ (no source 


correction), X-ray 


0.55 ±0.21 


7.23 


Lx, T x 


0.5 


X-ray only 




0.49 ±0.04 


3.40 


Y X ,T X 


2.5 


X-ray only 




2.19 ±0.16 


1.85 


Yx,Lx 


1.25 


X-ray only 




1.24 ±0.05 


4.29 



ables, rather than log-log space, and involved the minimi- 
sation of the squared and error-weighted residuals from a 
power law of the form y = Ax B , where A and B are to be 
estimated. The space of acceptable values of A and B was 
sampled by taking 10 5 datapoint realisations, distributed 
according to their errors, assumed Gaussian. The best-fit 
slopes are summarised in Table 6, with the values of the fit 
statistic. Figure 2(a) shows the scaling between yo and Tx. 
For this relation, the space of acceptable (A, B) parameters 
extends over a wide range in both parameters, so that the 
slope should be regarded as essentially unconstrained, al- 
though the best fit value B — 1.82 ± 0.72 is in agreement 
with the predicted value B = 1.5 from self-similar models. 
The other two relations show modest covariances between 
A and B and are more robust. Figure 2(b) shows the scal- 
ing between yo and Lx- We obtain a slope of 0.77 ± 0.21 
in good agreement with self-similarity, for which Equation 
6 predicts 0.75. Finally, Figure 2(c) depicts the scaling be- 
tween yo and Y%- We derive a slope of 0.65 ± 0.21 which 
is in good agreement with the predicted B — 0.6. For all 
three cases, the intrinsic scatter in yo can be estimated by 
the additional error required to bring the fit statistic down 
to a statistically acceptable level ~ 1.0 We find that the in- 
trinsic scatter in yo is about 25 per cent. The errors on B 
are large for all fits because of the modest range of cluster 
masses in our sample, but we see no evidence for departure 
from self-similar models of cluster evolution. Our results for 
the yo/Lx and yo/Tx scalings are consistent with those ob- 
tained by Morandi et al. (2007), although their gas masses 
and temperatures were defined over significantly different re- 
gions of the clusters (within -R2500 ) , and their scalings have 
smaller slope errors because of the wider range of Tx of their 
clusters. 

It is interesting to compare our scalings with those that 
we would have obtained if source corrections were not avail- 
able to us, since blank-sky surveys for cluster SZ effects 
could attempt to constrain the scalings with inadequate 
source data. We find that the scatter induced by the radio 
sources degrades the fit quality for each correlation signifi- 
cantly, as illustrated by the fit statistics shown in Table 6. 
Although the results broadly agree with self-similar predic- 
tions, the residual contamination tends to cause the scalings 



to appear too flat. For instance, in the high signal to noise 
sub-sample where the effect of sources is clearest, the slope 
of the yo/Lx relation is found to be flatter by 0.27 if a fit 
is attempted before source corrections are made. A similar 
result is found for the yo/Yx scaling relation. 

7.3 X-ray only scaling relations 

We check the consistency of the set of clusters with the more 
usual X-ray scaling relations between Lx and Tx, Yx and 
Tx, and Yx and Lx (e.g. Morandi et al. 2007) in the same 
way as for the SZ / X-ray scaling relations (results also 
summarised in Table 6). The slopes that we measure, of 
0.49 ±0.04, 2.19 ±0.16, and 1.24 ±0.05, are consistent with 
the similarity expectations of 0.5, 2.5 and 1.25. As the slight 
discrepancy for the correlation between Yx and Tx is not 
of high statistical significance, this study shows our sample 
to be representative of the population of hot (T > 6keV) 
clusters. 



8 DISCUSSION AND CONCLUSION 

We have observed a complete sample of galaxy clusters using 
OCRA-p, and studied the scaling of the central Compton 
parameter, yo, with various X-ray quantities. For each rela- 
tion, we find slopes in good agreement with the predictions 
from self-similar models. 

Our study has similarities to that of Morandi et al. 
(2007), and the two samples have 7 clusters in common, al- 
though we imposed stricter initial selection criteria. Morandi 
et al. consider scalings with yo, so we are able to perform a 
direct comparison. They find that the yo/Tx relation devi- 
ates by ~ 3a from the self-similar prediction, in the sense 
of being steeper. While our data are consistent with this re- 
sult, we note that they are essentially unable to constrain 
this relation at any level. Regarding the yo /Lx scaling, they 
find a slope of 0.61 ± 0.05, which is in good agreement with 
our work. Thus we see no sign of the flattening that would 
have resulted from the presence of an undetected set of con- 
taminating sources (Section 7.2). 

More recently, Bonamente et al. (2008) studied a sample 
of 38 clusters observed with OVRO/BIMA and Chandra. 
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They considered the integrated Comptonisation, Y, and in 
particular its scaling with Tx, M to t and M gaa . They find 
good consistency with predictions from self-similar models 
for all scaling relations. 

Our analysis has been limited by the narrow region of 
parameter space sampled, particularly for the scalings with 
Tx- Based on our results, a sample double the size would 
provide a good test of the yo/Vx scaling; this will be pre- 
sented in an upcoming paper in which we extend our sample 
to 33 clusters. To fully test the other relations, a sample four 
times larger is required. Thinking forwards to the prospects 
for blind cluster surveys, it is interesting to note that in this 
study we would have detected Sunyaev Zel'dovich effects for 
14 out of 17 clusters in our sample with no source informa- 
tion. Of course, to accurately determine the strength of the 
SZ effect in each case, and indeed to derive any cluster pa- 
rameters, good knowledge of the radio source environment 
is essential, and any attempt to test the scaling relations 
without taking sources into account will fail (see Section 
7.2). That this is true for sub-mm sources, as well as cm- 
wave studies, is clear from the large population of sub-mm 
sources known to be lensed by massive clusters (e.g. Johans- 
son et al. 2010). 
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